UTCCP-P-50, UTHEP-391 
September 1998 



1 



00 
On 

o\ 

Or 
00 



> 

OS 

o 

oo 

D 



X 



Quenched Light Hadron Spectrum with the Wilson Quark Action: 
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We report the final results of the CP-PACS calculation for the quenched light hadron spectrum with the Wilson 
quark action. Our data support the presence of quenched chiral singularities, and this motivates us to use mass 
formulae based on quenched chiral perturbation theory in order to extrapolate hadron masses to the physical 
point. Hadron masses and decay constants in the continuum limit show unambiguous systematic deviations from 
experiment. We also report the results for light quark masses. 



1. Introduction 

At Lattice'97 we presented first results from 
the CP-PACS calculation of the quenched light 
hadron spectrum with the Wilson quark action 
on large lattices (La>3 fm) at small quark masses 
(m^/mp = 0.75 down to 0.4) with high statis- 
tics (800, 600, 420 and 91 configurations at = 
5.9, 6.1, 6.25 and 6.47) 0. We have since in- 
creased the statistics at (3 = 6.47 to 150, and have 
completed the analysis. In this article, we report 
the final spectrum results and the main points of 
analyses behind them. 

2. Quenched chiral singularities 

Chiral extrapolation is a basic element of spec- 
trum calculations, for which a choice has to be 
made of the functional form of hadron masses 
in terms of quark masses. An important is- 
sue in considering the choice is the validity of 
quenched chiral perturbation theory(QxPT) ||- 
0, which predicts characteristic singularities in 
hadron masses in the chiral limit. We have there- 
fore made a detailed examination of this issue. 



2.1. pseudo-scalar mesons 

For pseudo-scalar (PS) mesons made of quarks 
of mass m and m s , QxPT formula reads |^] 



i 2 ps = A(m s + m){l - 5[\n{2mA/A 2 x ) 
+ m s /(m s — m) ln(m s /m)]} 
+ B(m s + m) 2 + C(m s — m) 2 + 



(1) 



To test the presence of the logarithm term, we 
combine our results to form two quantities 
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where it (rj) is the degenerate PS meson with 
quark mass m (m s ) at m^/rrip — 0.6, 0.5, 0.4 
(0.75, 0.7) and K is the non-degenerate one with 
m and m s . The two quantities are related by 
y = 1 + 5 ■ x, where the leading correction de- 
pends only on the 0((m s — m) 2 ) term in (Q). 

In Fig. we plot the two quantities calculated 
with quark masses determined from an extended 
axial current Ward identity (m^ WI ) as they have 
no ambiguity associated with the determination 
of the critical hopping parameter. The data fall 
within a narrow wedge spanned by the lines y ~ 
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Figure 1. Results of a test for the presence of 

the quenched chiral logarithms in PS meson mass 
data. See text for details. 

1 + (0.08 - 0.12)1, implying the value S w 0.10(2). 
We note that the 0((m s + to) 2 ) term can not be 
ignored for the range of our quark masses; results 
for the original ratio m^/mJ^J, which receive 
corrections both from 0{{m s +m) 2 ) and 0{{m s — 
to) 2 ) terms, do not fall on a common line. 

A different test using a ratio of decay constants 
y = /it /?)){§ leads to a similar result; our 
data fall within the lines y = 1 — 5/2 ■ x with 

5 = 0.08 - 0.16. 

Finally making full correlated fits to mpg us- 
ing ([!]) but imposing C = 0, independently for 
degenerate and non-degenerate data, we find 5 « 
0.06 - 0.12 for the range A x 0.6 - 1.4 GeV. 

These results lead us to conclude that our PS 
data show evidence for QxPT logarithms. 

2.2. vector meson and baryon masses 

For vector mesons and baryons, we perform 
uncorrelated simultaneous fits to degenerate and 
non-degenerate data together as a function of 
tops, assuming QxPT mass formulae |l||| with 

6 = 0.1. For vector mesons and decuplet baryons, 
all 0(mps) terms of QxPT arc included as well 
as 0(m PS ) terms. For octet baryons, we in- 
clude 0{mpg) terms in addition to 0(mps) and 
0(mp s ) terms since the nucleon mass shows a 
negative curvature which is opposite to that of 
the 0(mps) term. We omit decuplet-octet cou- 
pling (C in the notation of Ref. jfj)) and coupling 
to rj' (7), and set a$ = 0. 

These mass formulae fit our data well. The 
values for the coefficient C\i% of 0{mps) terms, 
however, are small. We obtain Cx/2 = —0.071(8) 
for p, -0.118(4) for nucleon, and -0.14(1) for A, 
to be compared with phenomenological estimates: 



1.8 
1.6 
1.4 



ffl 1.2 
CD 



1.0 
0.8 
0.6 
0.4 



K* 



1_ i 

T I Q 

1* 



□ GF1 1 K-input 
• CP-PACS K-input 
o CP-PACS $ input - 
experiment 



Figure 2. Final spectrum results compared to 
GFll's^l and experiment. 

C 1/2 = -AnglS « -0.71 for p if 8 = 0.1 and g 2 = 
0.75, and C 1/2 = -(3tt/2)(D - 3F) 2 S « -0.27 for 
nucleon if 6 = 0.1 and (F, D) = (0.5, 0.75). 

We conclude that vector and baryon masses are 
compatible with the presence of 0(mps) terms, 
but that their magnitudes are smaller than ex- 
pected. These terms yield an effect of at most 
0.2 x to^ w 25 MeV, which is about 3% of light 
hadron masses. 

3. Final spectrum results 

The results of analyses above motivate us to 
adopt the QxPT fits to calculate masses at each 
value of 13. The physical point for degenerate u 
and d quarks and the lattice scale are determined 
from the experimental values of to w and m p , and 
the strange quark mass by that of ttik or m^. We 
then extrapolate the results linearly in a. The fi- 
nal result for the spectrum in the continuum limit 
is shown in Fig j^. 

In order to examine how results differ if we 
do not employ Q%PT mass formulae, we repeat 
the analysis employing a quadratic polynomial 
in 1/K (cubic for N) for chiral extrapolations. 
While masses at each value of (3 differ, by about 
3% in the largest case, the differences in the con- 
tinuum limit do not exceed 1.5% of the results of 
QxPT fits. 

Compared to the results presented at Lat- 
tice'970] where we employed a linear chiral ex- 
trapolation in 1/K (cubic for N and quadratic 
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Figure 3. Light quark masses in MS scheme at 
[i=2 GeV. 

for A), the nucleon and A masses have decreased 
by 4.5% and 3.5%, respectively. Strange baryon 
masses with ttik used as input have also de- 
creased. The shift, however, is within 1.5<7 for 
all particles, with either ttik or as input. 

In summary, we find that differences in chiral 
extrapolations and an increase of statistics at (3 — 
6.47 do not alter the conclusions we drew at the 
time of Lattice'97: With m w , m p and ttik used as 
input, the meson hyperfine splitting and decuplet 
baryon mass splitting are too small compared to 
experiment, and so are the octet baryon masses. 
When we use instead of as input, the 
discrepancies for baryon masses are reduced, but 
the meson hyperfine splitting remains smaller. 

4. Light quark masses 

The Q%PT fit to pseudo-scalar meson masses 
has a significant effect on light quark masses at 
finite (3. Due to a negative curvature of the QxPT 
formula, values of the averaged u and d quark 
mass defined with vector Ward identity rn^Y 1 
become smaller than those from polynomial chiral 
extrapolations as shown in Fig. ||. The results 
extrapolated to the continuum limit, however, are 



consistent among various definitions. We adopt a 
combined fit to m^ WI and m^ WI , both estimated 
with QxPT fits, to calculate our final result. We 
obtain m U)d =4.6(2) MeV, and m a =115(2) MeV 
[m K input) or 143(6) MeV (m^ input) in the MS 
scheme at [i = 2 GeV. 

5. PS meson decay constants 

We determine f„ and Jk from the local axial 
current, employing a quadratic polynomial and 
linear chiral extrapolation, respectively. We ob- 
tain U = 120(6) MeV and f K = 139(4) MeV, 
which are 10 and 15% smaller than experiment, 
respectively. The ratio fx/ — 1 = 0.156(29) is 
also smaller than experiment. 

6. Conclusions 

We have presented our final results on the 
quenched light hadron spectrum and related 
quantities. In the course of analyses we found 
that our data for light hadron masses are con- 
sistent with predictions of QxPT. The effect of 
QxPT singularities is small, however, and the 
continuum results do not noticeably shift from 
those obtained with polynomial chiral extrapola- 
tions. Our results show that the quenched light 
hadron spectrum clearly and systematically devi- 
ates from the experimental spectrum. The dis- 
crepancy is much larger than our statistical error 
of 1% for mesons and 2-3% for baryons. 
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